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Abstract: We report the realization of a new iterative Fourier-transform 
algorithm for creating holograms that can diffract light into an arbitrary 
two-dimensional intensity profile. We show that the predicted intensity 
distributions are smooth with a fractional error from the target distribution 
at the percent level. We demonstrate that this new algorithm outperforms 
the most frequently used alternatives typically by one and two orders 
of magnitude in accuracy and roughness, respectively. The techniques 
described in this paper outline a path to creating arbitrary holographic atom 
traps in which the only remaining hurdle is physical implementation. 
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1. Introduction 



The intensity profile of optical fields can be controlled by using a computer-generated digi- 
tal phase raster, also called a computer-generated hologram (CGH). The raster, or kinoform, 
changes the complex phase of the optical field at each pixel without attenuating the amplitude. 
The CGH, physically implemented using a spatial light modulator (SLM) or similar technology, 
is illuminated by a monochromatic beam which is relayed to a focusing objective; the desired 
intensity distribution is produced at the focal plane of the objective. CGHs have been used in 
a wide range of technological applications and physical research. For example, soft matter and 
biological systems have been manipulated using dynamic arrays of focused beams (i.e., optical 
tweezers) created using CGHs (see 1 1 1 for a review), and CGHs have been used to shape laser 
beams for inertial confinement fusion experiments |l2l[3]|4l. 

Although CGH technology has recently been applied to experiments involving ultra-cold 
atom gases, that work has been limited to arrays of optical dipole traps fSl |6l H [H. New ap- 
plications of optical dipole potentials for cold-atom research, such as interferometers ||9l and 
atom "transistors" |10|, may be realized if high-quality arbitrary light intensity profiles could 
be generated. The trapping potential for an atom confined in an optical dipole trap is propor- 
tional to the light intensity, and therefore the properties of the intensity profile created by a 
CGH are critical for such applications. Desirable features of a CGH applied to trapping ultra- 
cold atoms include accuracy in matching the desired intensity profile, smoothness of the profile 
generated by the CGH, and efficiency in diffracting light into the target profile. Smoothness, of 
particular importance for ultra-cold atom experiments, has received limited attention in previ- 
ous work on CGHs lflTl[T2l . Disorder related to small-length-scale intensity fluctuations in an 
optical potential will introduce complications for interferometric applications lfT3l [141 [T5l [T6l 
and can greatly affect the study of quantum many-body physics (see ifTTl . for example). Pre- 
vious measures used to analyze the roughness of intensity profiles created using CGHs could 
only be applied to uniform distributions; in this manuscript, we introduce a new roughness met- 
ric appropriate for continuous profiles. Some constraints related to CGHs may be relaxed for 
ultra-cold atom applications — a CGH in this context is not required to control the intensity in 
the entire focal plane, since ultra-cold atom gases are typically confined to a finite region of 
space. 

Calculating a CGH to generate a high-quality arbitrary light intensity distribution is a chal- 
lenging problem, because a CGH cannot be directly computed, in general, from a desired arbi- 
trary intensity profile. One technique for calculating a CGH when an exact solution is unknown 
is to use an iterative Fourier transform algorithm (IFTA), which is computationally efficient 
compared with other methods, such as a direct binary search ifTSl [T9l . An IFTA predicts the 
propagation of a beam through an initial kinoform by fast Fourier transform (FFT), and then 
successively modifies the kinoform based on a comparison between the predicted and desired 
focal plane intensities. The most frequently used IFTAs for calculating CGHs are variants of 
the Gerchburg-Saxton (GS) and Adaptive-Additive (AA) algorithms fm l20l IITI l22l l23l . 

In this paper, we present a new IFTA that we call the "mixed-region amplitude freedom" 
(MRAF) algorithm. The MRAF algorithm typically improves by one order of magnitude on 
accuracy and and two orders of magnitude on roughness compared with the GS and AA algo- 
rithms for continuous target profiles. To our knowledge, no algorithm for creating CGHs with a 
comparable level of computational complexity surpasses the MRAF algorithm in measures of 
accuracy and roughness. The MRAF algorithm controls intensity in a bounded two-dimensional 
subset of the focal plane and achieves accuracy at the percent level at typically the cost of a 
factor of three in efficiency (compared with the GS and AA algorithms). Because the MRAF 
algorithm controls the intensity profile in a single plane, this method can only be applied to 
creating two-dimensional arbitrary optical traps; confinement to the focal plane will require an 




Fig. 1. Schematic geometry for an IFTA. Tlie optical field that propagates from the input 
to the output plane through a focusing objective is shown in green. The field is discretized 
using coordinates {x,y) in the input plane and {x' ,y') in the output plane. The dashed lines 
represent the clear aperture of the focusing optics. The matrix used to computationally 
represent the input field must be enlarged beyond this region and filled with zero intensity 
points (dark gray) to fully resolve the output plane. The physical size of the matrix used to 
represent the input plane is d. 



additional tightly-focused sheet of light. In section [2] of this manuscript the MRAF algorithm 
is explained in detail, and in section |3] we report on the algorithm performance for six target 
intensity profiles. 

2. Algorithm 

Before giving the mathematical details of the MRAF algorithm, we briefly review the operation 
of an IFTA [24]. An IFTA is a technique to solve the following problem: design a CGH that 
will convert a light field AQ{x,y) at the CGH, or input plane, into a target intensity distribution 
/o(x',y) at the focal, or output, plane of a focusing optic (see Fig.[T]i|25|. The light field Aq is 
typically a Gaussian beam apodized by the input aperture of the device used to implement the 
CGH. The IFTA problem does not have a unique solution, as the complex phase of the optical 
field associated with Iq is not constrained. This is known as phase freedom; there is a choice 
of phase in the output plane 1261 1271 . Complete phase freedom is allowed for far-off resonance 
optical atom traps because the phase of the field in the output plane does not contribute to 
the trapping potential if the light is far-detuned from an electronic transition and if the dipole 
approximation is valid. An IFTA is designed to use phase freedom to minimize the difference 
between Iq and the intensity distribution produced by the CGH in the output plane. 

An IFTA can be decomposed into two parts as illustrated in Fig. [2] an initialization step and 
an iterative loop. In the initialization step, a phase distribution Ko{x,y) is chosen as a starting 
point for the algorithm and is imprinted on Aq to produce the input field E^^\x,y) = |Ao|e'*^ 
for the first iteration. Each iteration n of the loop begins by calculating the field Ei",l{x',y') = 
^ Ej"^ produced by E^""^ propagating to the output plane. The propagation is modeled using 
a Fourier transform which assumes the paraxial approximation for the focusing optics ||28ll . 

(n) 

The algorithm then combines the propagated field E^J with the target intensity profile Iq to pro- 
duce a new field G''"'(x',y). This procedure is carried out using one or more numerical scalars 
called mixing parameters m. The phase of the backward propagated field arg [^^' [G(n)]] is 
used as the starting phase distribution for the next iteration. 

The iterative loop is terminated after iterations once a figure-of-merit 7], calculated using 
the intensity profile in the output plane and Iq, does not improve with repeated iterations. The 
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Fig. 2. Block diagram of an IFTA. 



phase profile Kf{x\y') = arg of the field in the input plane for the final iteration is 

the kinoform which must be transferred to a physical device. An ideally implemented CGH 
will produce the predicted intensity profile If{x',y') — \eI,^)\^ in the output plane. The goal 
of using an IFTA to design a CGH is to choose mixing parameters that optimize one or more 
measures calculated on the predicted profile //, such as the deviation from Iq. Achieving this 
goal will typically require executing an IFTA multiple times with different selections of mixing 
parameters as part of an optimization scheme |29|. 

Central to the MRAF algorithm is the introduction of amplitude freedom into a restricted 
region of the output plane and the use of a single mixing parameter. The mathematical details 
of the MRAF algorithm are given in section 2.1 The use of amplitude freedom in the MRAF 
algorithm is not sufficient to generate a high-quality optical field — a choice of initial phase that 
eliminates optical vortices from the output plane for all iterations is also necessary. We outline 



a procedure for choosing appropriate initial phase distributions in section 2.2 



2.1. MRAF algorithm 

At each step n of the MRAF algorithm, the propagated field is combined with the target intensity 
distribution according to: 

gW ^ {m^^sR + {l-m) \E^\Je'<'-'\. (1) 

A single mixing parameter m controls the relative distribution of optical power in two subsets, 
the signal region {SR) and noise region (A^^), of the output plane. Phase freedom is permitted 
everywhere in the output plane (the phase of the propagated field is used for the phase of 
G'"^), while amplitude freedom is allowed only in the noise region. Even though the mixing 
parameter m is kept fixed, the fraction of power in each region changes for every iteration; the 
only constraints imposed on the MRAF algorithm are that the power in the target profile ^/q 



(which is only non-zero in the SR) and the total power in the output plane niY^sR 
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remam constant. 

The signal region is chosen to overlap with the area in which light will interact with atoms; 
the remainder of the output plane is the noise region. The effect of dividing the output plane into 
subsets is to cause the algorithm to converge very closely to the target profile within the signal 
region, while behaving in a less controlled manner in the noise region. Utilizing the amplitude 
freedom in the noise region allows for increased accuracy in matching Iq in the signal region, 
while decreasing the efficiency of the CGH. The MRAF algorithm is equivalent to a variable 
strength projection II30I or regularized algorithm 131] with a specific trajectory for the variable 
mixing parameter that has not been previously demonstrated. The MRAF algorithm is also 
similar to the algorithm used in |27|, but with a fixed m. 

In Sec.|3] we compare the performance of the MRAF algorithm with the GS and AA algo- 

rithms. The GS algorithm, in which G'"^ — ^/loe'^^i °"'J, permits only phase freedom in the 
output plane. In the AA algorithm, amplitude freedom is introduced uniformly into the output 

plane: G^ = {m^+ (1 - m)£i"j} . 
2.2. Initial phase 

As in any optimization scheme, an initial guess that produces a result close to the target im- 
proves the convergence rate and reduces the risk of stagnation into a local optimization min- 
imum. Because IFTAs are used when the solution to the CGH problem is unknown, choos- 
ing an initial phase profile Kq to reproduce complex features in Iq is not possible. Therefore, 
we wish to find a Kq as our starting point for which most of the power in the output plane 
roughly overlaps with the envelope of Iq. The distribution Kq must also be chosen so that 
does not contain any undesired optical vortices — points characterized by a phase singularity 
and zero intensity — since an IFTA is not able to eliminate vortices present in the output plane 
ll32l [33l l34l . A further constraint is that only certain choices for Kq can prevent an IFTA from 
producing optical vortices in the output plane at each iteration 132] [33l [34l . The optical field 
in an IFTA is discretized and the Fourier transforms are calculated using FFTs. In order for 
the FFT to fully resolve the output plane and to reproduce the aperture of the physical device 
used to implement the CGH, the matrix representing the optical field must be enlarged (i.e., 
"padded") beyond the size of the CGH by adding a zero-intensity region lf35l . The truncation 
of the input plane field caused by padding leads IFTAs to create optical vortices in the output 
plane; this behavior is not completely understood [32] . 

We find that a quadratic phase distribution for Kq combined with linear and conical gradients 
does not introduce optical vortices for the MRAF algorithm, even through the input field is 



truncated in each iteration of the IFTA. Quadratic phase distributions, the equivalent of a thin 
lens, were first discussed as a solution to the vortex problem in the context of variants of the GS 
algorithm 1321 . A quadratic phase profile, given by Ko{x,y) = 4-R \ax^ + (1 — Ct)}'^] (R is the 
curvature and a/(l — a) is the aspect ratio) changes the size of the envelope of the intensity 
profile in the output plane. A linear gradient phase profile K{){x,y) = B [xcQ^,{pL) +ysm{pL)], 
where B is the strength of the gradient and pL is an angle, shifts the centroid of the intensity 
profile in the output plane. A conical phase gradient Koix^y) — Br creates a ring in the output 
plane, where r = \J x^ +3'^. 

For the work in this paper, quadratic phase profiles are used to roughly match the size of the 
field in the output plane to the size of the target profile. We find that the results of the algorithm 
are not strongly affected by small changes in R. Linear gradients are used to match targets that 
are shifted from the center of the output plane to avoid complications caused by undiffracted 
light resulting from the finite efficiency of a physically-implemented CGH. Conical gradients 
are used to match target profiles which have a ring-Uke structure. Quadratic, linear, and conical 
phase distributions are added together modulo 2n to combine the effects of each. We find that 
this combination of initial phase distributions has enough flexibility to obtain sufficient overlap 
of \E^out P with /q to achieve a few percent error in If within tens of iterations. 

3. Results 

We characterize the performance of the MRAF algorithm using the six target intensity profiles 
shown in Fig. |3] chosen because of their potential appUcation to ultra-cold atom experiments. 
Target (a) consists of two Gaussian beams connected by a ring with a Gaussian cross-section, 
which could be used to study ultra-cold atom gases in multiply-connected geometries. A 3- 
pointed star-shaped intensity profile is shown in Fig.jSj^b); a similar profile was recently used 
to induce spontaneous vortex generation in an atomic Bose-Einstein condensate (BEC) lf 36ll . 
Fig. |3jc) shows a uniform square intensity profile, which may be used as an optical lattice 
beam in experiments for which it is desirable to remove effects generated by the external con- 
finement resulting from a Gaussian beam profile f37' '38l. Target (d) is a complex intensity 
profile designed to be evocative of an optical lattice beam that could be employed to realize an 
"atomtronic" logical OR gate ifTOl . An intensity profile that could be used to trap a BEC in a 
superconducting quantum interference device (SQUID) geometry is shown in Fig.[3je). Finally, 
a BEC confined in a dipole trap created using target profile (f) would be equivalent to a thin 
superconducting wire connected between bulk superconductors |f39l . 

The parameters defining the geometry of each target profile are given in the caption to Fig. [3] 
To describe these test patterns and in the rest of this manuscript, we use pixels (px) to measure 
distances in the input and output planes. Each pixel in the input plane represents a point at which 
the input field is discretized. The physical size of a pixel in the input plane is the physical size, 
d, of the matrix used to represent the input field divided by the number of pixels. The pixel size 
in the output plane is /A jd, where / is the focal length of the focusing objective and X is the 
wavelength of the light. For the results given in this section, the input field is discretized on a 
768 X 768 pixel array and the phase of the field in the input plane is discretized in 256 levels 
(for each iteration of the IFTA) ll40l . The input field array is padded with zero intensity points 
in each iteration of the IFTA to create a 1536 x 1536 pixel matrix. We observe no significant 
change in the algorithm results if the array is enlarged beyond 1536 x 1536 pixels (consistent 
with the Nyquist-Shannon sampling theorem). We implement IFTAs in MATLAB 141 1, and we 
use a Gaussian input field Aq °= e^' ' "0 with a 565 pixel waist wq. 

The initial phase profile chosen for each target intensity profile is shown in the top half of 
Fig. |4] A conical phase profile is used in for (a) and (e) to create a ring structure. Target 
profiles (b) and (c) are shifted from the center of the output field, and therefore linear gradients 



(a) 

< 1 


(b) 

Y 


(c) 

■ 


(d) _ 

1- 


(e) 1 




• 



Fig. 3. Target intensity profiles Iq used to ciiaracterize the performance of the MRAF algo- 
rithm. The field-of-view for images (a), (b), (c), and (e) is a 200 x 200 pixel and for (d) and 
(f) is a 400 X 400 pixel subset of the output plane centered on /q. The grayscale represents 
intensity, with black corresponding to the regions of zero intensity. The radius of the ring in 
(a) is 53 pixels and the waist for each Gaussian beam and the Gaussian cross-section of the 
ring is 14 pixels. The maximum intensity of the Gaussian beams is three times that of the 
ring. Each "tip" of the star-shaped pattern in (b) is 20 pixels from the center of the star; the 
two lines that intersect to form each "tip" subtend a 28° angle. To create the profile shown 
in (b), a uniform intensity profile with a star shape was convolved with a Gaussian with a 5 
pixel waist. The profile in (c) was created by convolving a 58 pixel on edge square profile 
with a 3 pixel averaging filter. The overall dimensions of the profile in (d) are 288 pixels 
wide and 325 pixels high, and the intensity in the "base" regions is increased by 33%. The 
Gaussian ring in (e) has a 53 pixel radius and a 7 pixel r.m.s. width. The intensity in the 10 
pixel wide gaps in (e) is suppressed by a factor of 2, and the "leads" in (e) are 185 pixels 
from end-to-end. The 264 pixel wide Gaussian "wire" in (f) has a 3.5 pixel r.m.s. width, and 
the Gaussian reservoirs in (f) have a 17.6 pixels r.m.s. radius. The center of each profile is 
displaced from the center of the output plane by (b) 37 pixels and (c) 63 pixels; the profiles 
in (a), (d), (e), and (f) are centered on the output plane. 



are used in Kq to displace the intensity in ii^^'j accordingly. Quadratic phase profiles are used in 
each Kq to match the approximate area covered by Iq in the output plane. The intensity profile 
in the output plane for the first iteration of the IFTA is shown in the bottom half of Fig.|4] The 
initial phase profiles in Fig. |4] were optimized manually; the predicted intensity profile is not 
affected by small changes in A'q. 

A qualitative comparison between the predicted intensity profiles for the MRAF, GS, and 
AA algorithms is shown in Fig. |5] We find that the small-length-scale intensity fluctuations 
apparent in Fig.|5]are a generic feature of using the GS and AA algorithms to generate arbitrary 
intensity profiles. We did not determine if any of these fluctuations are optical vortices, which 
can be removed from If under limited circumstances by changing Kf (and at the cost of greatly 
increased computational complexity) ll22l |32l |34l |42 | . 

To quantitatively characterize the performance of the MRAF algorithm and compare with 
the GS and AA algorithms, we introduce measures of accuracy and roughness. The accuracy 
metric is the root-mean-square (rm.s.) fractional error from the target averaged across a subset, 
the measure region (MR), of the output plane: 



' E 



ifix',y)-iQ{x',y) 
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The measure region is a subset of the signal region and is chosen to exclude the zero-intensity 
pixels in Iq. The intensity profiles // = ^f/'L{x'.y')eMR^f ^^d Iq — Io/'L{x',y')eMR^o normal- 
ized to have the same power in the measure region, and Nmr — Jl{x'.y)eMR 1 is the number of 
pixels in the measure region. The error rj is also used as the convergence parameter and the 
optimization parameter for the MRAF and AA algorithms. 
The roughness measure 

P= E {H\rfix'y)-k{x',y)]y /Nmr (3) 
{x'y)eMR 

is the average of the square of the mean curvature H of the difference between the predicted 
and target intensity profiles in the measure region. The roughness p is proportional to the Will- 
more bending energy for the surface Iq — If |43 1. A unique measure of roughness for a two- 
dimensional manifold, such as //, does not exist. We choose p as defined in Eq. 3] as a metric 
because it is intuitively appealing as an energy that is strongly weighted by smal -length scale 
deviations of If from Iq (the Willmore bending energy of a spherical surface is proportional 
to the square of the inverse of the radius of curvature). The measure p also quantitatively re- 
produces qualitative features we observe in the predicted intensity profiles. For each predicted 
intensity profile we also calculate the efficiency of the CGH for diffracting light into signal re- 
gion. The efficiency ^ — L(T',y)G5/j^//L(A',>')^/ is defined as the ratio of the power in the signal 
region to the total power in the output plane. The parameters rj and p are both efficiency- 
independent measures of the deviation of If from Iq: rj and p are zero if If = Iq. 

The result of using the MRAF algorithm to create a CGH based on the test target profiles 
is shown in Fig. |6] Both the final kinoforms to which the the initial phase profiles in Fig. |4] 
converged and the predicted intensity profiles are shown in the figure. The mixing parameters 
and the conical and quadratic phase profiles used in Kq were optimized by determining the 
minimum value for rj calculated for a wide range of m, B, and R. The signal region used in 
the MRAF algorithm is outlined in red for each target profile in Fig.|6] The MRAF algorithm 
converged in less than 100 iterations for each of these target profiles. 

Table 1 shows the accuracy, roughness, and efficiency calculated for each predicted inten- 
sity profiles in Fig. |6] and the equivalent results generated by the GS and A A algorithms. The 
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Fig. 4. Initial phase profiles Kq (top row) and predicted initial intensity profiles l-Eoi) P (bot- 
tom row) chosen for the target intensity profiles in Fig.|3] The phase profiles are 768x768 
pixels and are shown in grayscale modulo 2n; white corresponds to a In phase. Conical 
phase profiles with 6 = 117 mrad/px are employed in (a) and (e). Linear gradients of 
136 mrad/px and 260 mradj px with /i = and n/A are used for (b) and (c), respectively. 
Quadratic phase profiles withi? = 0.31 mrad/px^, R = 0.3 mrad/px^, if = 0.34 mrad/px^, 
i? = 1.4 mrad/px^, R = 0.5 mrad/px^, and R= 1.6 mrad/ px^;a = 0.29 are applied in (a), 
(b),(c), (d), (e), and (f) respectively. 
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Fig. 5. Intensity profiles If for the MRAF, GS, and AA (left, middle, right) algorithms for 
the test target profiles. Only the intensity in the signal region is shown, and the profiles are 
scaled so that the total power in the signal region is the same for each. 




Fig. 6. Final kinoforms Kf (top row) and predicted intensity profiles If (bottom row) pro- 
duced by the MRAF algorithm for targets (a), (b), and (c). The mixing parameters used to 
generate these results are: (a) 0.40, (b) 0.35, (c) 0.40, (d) 0.30, (e) 0.35, and (f) 0.30. The 
signal region (red) for (a) is an annulus with with inner and outer radii 25 and 81 pixels; 
in (b) is a circle with a 40 pixel radius; for (c) is a square 75 pixels on edge; in (d) is a 
region 10 pixels from the edge of the target profile; for (e) a region 10 pixels away from 
where the intensity is 10% of the maximum intensity; and for (f) consists of two 53 pixel 
radius circles separated by 264 pixels and a connective region 25 pixels wide. The measure 
region (yellow) in (a) is an annulus with inner and outer radii 44 and 62 pixels; in (b), (e), 
and (f) is defined by a region in which the intensity of the target is greater than 10% of the 
maximum target intensity; in (c) is a square 57 pixels on edge; and for (d) is the edge of the 
target profile. 



accuracy and roughness are determined using the measure region outlined in yellow in Fig. [6] 
The MRAF algorithm on average shows a factor of 9 improvement in accuracy and a factor of 
190 improvement in roughness compared with the GS and AA algorithms; the average error 
for each target for the MRAF algorithm is at the few percent level. The MRAF algorithm pro- 
duces comparatively smooth intensity profiles even though p is not used to optimize the mixing 
parameters or as a convergence criteria for the IFTA (see Fig. [7]). 



Algorithm 


Error rj 




Ring (a) 


Star (b) 


Square (c) 


OR gate (d) 


SQUID (e) 


thin wire (f) 


GS 


0.21 


0.30 


0.23 


0.34 


0.36 


0.36 


AA 


0.13 


0.19 


0.19 


0.23 


0.21 


0.23 


MRAF 


0.017 


0.027 


0.015 


0.039 


0.018 


0.029 




Roughness p 


GS 


220 


5600 


460 


13 


160 


110 


AA 


64 


2200 


400 


5.3 


47 


40 


MRAF 


0.65 


20 


1.1 


0.044 


0.18 


0.19 




Efficiency ^ 


GS 


0.99 


0.99 


0.99 


0.96 


0.97 


0.97 


AA 


0.90 


0.83 


0.88 


0.79 


0.71 


0.59 


MRAF 


0.45 


0.29 


0.45 


0.18 


0.30 


0.19 



Table 1. Table comparing the performance of the MRAF to the GS and AA algorithms. The 
mixing parameters used for the AA algorithm are (a) 1.9, (b) 2.0, (c) 1.9, (d) 2.0, (e) 2.2, 
and (f) 2.5. The GS and AA algorithms converged in 100 iterations for the results in this 
table. 

The error Tj, an average across the measure region, can be small although large errors exist at 
points in the output plane. To show that the MRAF algorithm achieves accuracy everywhere in 

the signal region, a histogram of the error evaluated at each output plane point {\J l/' ~Iq^ /Iq) 
in the signal region for the result in Fig.|6|a) is shown in Fig. [8} 95% of the pixels in the signal 
region have less than a 3% error. For the purposes of comparison, the result of using the GS and 
AA algorithms to calculate a CGH for target (a) is also shown in Fig.|6] The MRAF algorithm 
improves greatly on the GS and AA algorithms, both of which have at least 45% of the pixels 
in the signal region with errors greater than 10%. 

For the targets used in this manuscript, the MRAF algorithm has approximately a factor of 
2-3 lower efficiency compared with the GS and AA algorithms. While the MRAF algorithm 
does not lead to efficiencies as high as the GS and AA algorithms, it does not perform so poorly 
that the efficiency advantage of using a MRAF is lost. For example, the authors of |36| created 
a profile similar to (b) by propagating a Gaussian beam through an intensity mask and imaging 
the result onto a BEC. The efficiency of that mask in transmitting light was 3%, compared with 
« 29% for a CGH created using the MRAF algorithm. 

4. Conclusion 

In conclusion, we have reported the realization of a new IFTA for designing CGHs that can 
be used to create arbitrarily-shaped, two-dimensional optical dipole traps for ultra-cold atom 
experiments. The MRAF algorithm has relatively low computational complexity and converges 
rapidly — within tens of iterations. For six test target profiles, the predicted output of a CGH 
designed using the MRAF algorithm is comparatively smooth and has errors at the percent 
level. 
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Fig. 7. Variation of measures characterizing the MRAF algorithm performance as the 
mixing parameter m is varied. The efficiency ^ (red), roughness p (black), and error 77 
(blue) are shown for target (a) for different values of the mixing parameter m. The inset 
shows detail around the globally-optimized value of m. The mixing parameter that mini- 
mizes rj approximately coincides with an minimum in roughness p for the MRAF algo- 
rithm. 

The challenge for realizing arbitrary, two-dimensional dipole traps for atoms now lies with 
experimentally implementing a CGH designed using the MRAF algorithm. In this paper, we 
assume abberation-free optics, the paraxial approximation, a single-polarization optical field, 
and ideal CGH response. Some or all of these idealizations will be violated in an experimen- 
tal realization of a CGH, leading to the intensity profile produced at the focal plane deviating 
from the predicted profile If. The extent to which these practical considerations affect an ex- 
perimental implementation will depend on the details of the kinoform and the specific nature 
of technical problems. 

In particular, non-ideal CGH response is likely to have a high-impact on the quality of the 
intensity profile produced in an experiment. In this paper, we take into account two practical 
limitations of CGH technology: quantized phase levels and finite resolution. Many approaches 
to producing CGHs also do not have well-characterized or well-controlled phase response. 
For example, commercially-available, scientific-grade SLMs are afflicted by non-uniform and 
nonlinear phase response |44| . Controlling these problems at the percent level to experimentally 
achieve the high-accuracy of the MRAF algorithm will be challenging. 
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Fig. 8. Histogram of the fractional error at each pixel evaluated for If for the MRAF algo- 
rithm used on target (a). The fraction of pixels in the signal region are binned with respect 

/~2 ~2 ~ 

to the fractional error y If —Iq /Iq. The width of each bin is equivalent to a 1% fractional 
error. The solid black, blue dotted, and red dashed lines are the result for the MRAF, GS, 
and AA algorithms, respectively. 



